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(N : ABSTRACT 

We develop an efRcient method to study the effects of reionization history on 
> ■ 

' the temperature-density relation of the intergalactic medium in the low den- 

m ■ 
(N 



sity limit (overdensity 6^5). It is applied to the study of photo-reionization 
models in which the amphtude, spectrum and onset epoch of the ionizing flux, 
i as well as the cosmology, are systematically varied. We find that the mean 

o\ 

' temperature-density relation at z = 2 — 4 is well approximated by a power- 

Oh' law equation of state for uniform reionization models. We derive analytical 

O I expressions for its evolution and exhibit its asymptotic behavior: it is found 

C/3 . that for sufficiently early reionization, imprints of reionization history prior 

to 2; ~ 10 on the temperature-density relation are washed out. In this limit 
' the temperature at cosmic mean density is proportional to [Qbh/y/QQ]^^^'"^ . 

^ I While the amplitude of the radiation flux at the ionizing frequency of H i is 

found to have a negligible effect on the temperature-density relation as long 
as the universe reionizes before z ~ 5, the spectrum can change the overall 
temperature by about 20%, through variations in the abundances of helium 
species. However the slope of the mean equation of state is found to lie within 
a narrow range for all reionization models we study, where reionization takes 
place before z ~ 5. We discuss the implications of these findings for the ob- 
servational properties of the Lya forest. In particular, uncertainties in the 
temperature of the intergalactic medium, due to the uncertain reionization 
history of our universe, introduces a 30% scaling in the amplitude of the col- 
umn density distribution while the the slope of the distribution is only affected 
by about 5%. Finally, we discuss how a fiuctuating ionizing field affects the 
above results. We argue that under certain conditions, the loss of memory 
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of reionization history implies that at late times, the temperature-density re- 
lation of a gas in a fluctuating ionizing background can be approximated by 
one that results from a uniform radiation field, provided the universe reionizes 
sufhciently early. 

Key words: cosmology: theory — intergalactic medium — quasars: absorp- 
tion lines 

1 INTRODUCTION 

Absorption line studies of the low column density Lya forest (A^hi ~ 15cm~^; see Hu et al. 
1995; Lu et al. 1996; Cristiani et al. 1996) offer a unique probe of the universe which is 
free of luminosity bias, spans a wide redshift range z ~ 2 — 4 and where the overdensity is 
mildly nonlinear (overdensity 5^5, see Hui, Gnedin & Zhang 1996) and hence the relevant 
physics is simple to understand. Recent hydrodynamic simulations support the hypothesis 
that the Lya forest arises from a fluctuating intergalactic medium which is the natural result 
of gravitational instability in hierarchical clustering cosmological models (Gen et al. 1994; 
Zhang et al. 1995; Hernquist et al. 1996; Miralda-Escude et al. 1996; Wadsley & Bond 1996). 

The main quantity of interest in such studies is the Lya transmission (or optical depth) as 
a function of frequency, which is determined by the density of neutral hydrogen as a function 
of position. The Gunn-Peterson effect (Gunn & Peterson 1965) tells us that hydrogen (as 
well as helium) is highly ionized for z^5. Several models for its origin have been studied: 
photoionization by first generation of stars or quasars, collisional ionization induced by 
shock-heated gas, decaying neutrinos and so on (Ikeuchi & Ostriker 1986; Gouchman & 
Rees 1986; Sciama 1990; Miralda-Escude & Ostriker 1992; Shapiro, Giroux, & Babul 1994; 
Giroux & Shapiro 1996; Haiman & Loeb 1996). The reader is also referred to Shapiro (1995) 
for an excellent review on the subject. We focus on photoionization models in this paper 
although the semi-analytical method we propose in this paper can be applied to any model. 

Ionization equilibrium, which should hold after the universe reionizes, dictates that the 
neutral hydrogen density is approximately proportional to pfT~°'^ where pb is the baryon 
density and T is the temperature (see § ^ . 1 . 1| ) . It is therefore important to know the am- 
plitude of T and its dependence on pb. As we will discuss later, the temperature-density 
relation can significantly affect observed properties of the Lya forest, its column density 
distribution for instance. 

Imagine a fluid element evolving in a photoionized intergalactic medium. Its temperature 
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at ^ ~ 2 — 4 is influenced by a number of factors: how its density evolves with time, which 
determines the amount of heating/coohng; how and when the ionizing radiation is turned on; 
the radiation amphtude and spectrum, which control how much photoionization heating the 
fluid element suffers; cosmological parameters such as baryon density, the Hubble constant 
and the density parameter which affects the recombination rate, adiabatic cooling rate 
and the evolution of density. 

The temperature-density relation for a set of fluid elements is then determined by, among 
other things, the reionization history of the universe. Given a reionization scenario and cos- 
mological model, hydrodynamic simulations can predict this relation accurately, but limited 
computer resources obviously restrict the number of reionization histories one can study. It is 
therefore important that we develop alternative tools in order to understand the effects of a 
full range of reionization histories, considering the as yet poor knowledge of the actual reion- 
ization history of our universe. Moreover, it is worth pointing out that there exists analytical 
tools borrowed from studies of large scale structure which can predict the distribution of pb 
with reasonable accuracy in the low density regime (Bi, Borner & Chu 1992; Bi & Davidsen 
1996; Gnedin & Hui 1996; Hui, Gnedin & Zhang 1996). By their analytical/semi- analytical 
nature, it is possible to study a whole range of cosmological models in an efficient manner. 
What remains to be specified for these methods, in order to calculate quantities hke the 
column density distribution of hydrogen, is the temperature-density relation. A calculation 
aiming at just that is attempted here. 

We propose to model the density evolution of each fiuid element using the Zel'dovich 
approximation (Zel'dovich 1970), which is known to be a good approximation in the mildly 
nonlinear regime (Coles, Melott, & Shandarin 1993) and for elements where hydrodynamic 
effects like pressure and shock- heating are unimportant. Given a probability distribution of 
the initial configuration for these elements, which can be deduced for any given cosmological 
model, one can generate a whole set of such elements and simply follow their thermal and 
chemical evolution. The resulting bulk temperature-density relation can then be studied 
for any chosen reionization history. The virtue of this method is that each element can be 
treated independently and that the density evolution is given by an analytical formula (from 
the Zel'dovich approximation) while the thermal and chemical evolution equations involve 
a set of ordinary differential equations which can be solved numerically. Obviously, such a 
Lagrangian method does not incorporate the effects of shock-heating properly. We check the 
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results of our method against full hydrodynamic simulations and find good agreement for 
5^5, implying shock effects are not significant in the low density regime (§ |2.2|) . 

Most previous work investigating the effects of different reionization histories either as- 
sumes a uniform intergalactic medium (eg. Giroux & Shapiro 1996) or uses the spherical 
collapse model for density evolution (eg. Miralda-Escude & Rees 1994). For the low column 
density Lya forest, which arises naturally from a mildly fluctuating intergalactic medium, 
and which hydrodynamic simulations indicate to be consisting mostly of filaments and pan- 
cakes (Gen et al. 1994), the density evolution is better approximated by the Zel'dovich 
approximation (Hui et al. 1996). 

A great deal of effort has been made to examine radiative transfer effects after the onset of 
ionizing radiation, including the expansion of H ll regions, reprocessing by absorption systems 
etc. (Zuo 1992a,b; Meiksin & Madau 1993; Miralda-Escude & Rees 1994; Giroux & Shapiro 
1996; Haardt & Madau 1996; see also Gnedin & Ostriker 1996 for numerical simulations that 
include some of these effects). These effects can in principle be incorporated in our method 
by specifying how the radiation intensity correlates with density and possibly allowing the 
onset of reionization to occur at different times for different fluid elements. While more work 
still needs to be done in order to make such a speciflcation with reasonable accuracy, we 
focus in this paper on reionization models in which the ionizing flux is uniform, but allowing 
the onset-epoch, spectrum and amplitude of the ionizing flux as well as the cosmology to 
vary. It turns out one can learn something about the more realistic fluctuating case, even 
from these simple models, which we will discuss in § |[ 

As we will show, the mean temperature-density relation of the photoionized intergalactic 
medium in the low density regime is well-approximated by a power law equation of state 
(T oc [1 + 5]'''^"'^) where 6 is the mass overdensity) . The questions we would like to address are 
then: how does the amplitude and slope of the equation of state depend on the reionization 
history, characteristics of the ionizing radiation as well as cosmological parameters; what are 
the implications for observational properties of the Lya forest? 

Organization of the paper is as follows. In the next section, we outline a semi-analytical 
method to predict the temperature-density relation for any given cosmology and reionization 
history, using the Zel'dovich approximation. An example is discussed and a comparison is 
made with the result of a full hydrodynamic simulation, computed using the same reion- 
ization history. Then, in § |] we apply the method to study models in which the universe 
reionizes suddenly, due to a rapid outburst of radiation. We investigate systematically how 
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varying the epoch of reionization, cosmological parameters and the radiation amphtude and 
spectrum affects the equation of state at later redshifts. Analytical approximations to the 
equation of state are derived. In § ^, it is shown that models in which reionization occurs 
after a period of reheating (gradual turn-on of radiation) give equations of state that can 
be well-fitted by that of an appropriate sudden reionization model. Finally we discuss the 
implications of our findings for the observational properties of the Lyman-alpha Forest and 
the intergalactic medium in § ^. A summary of relevant reaction rates is provided in the 
Appendix. 

A few words about our terminology. The three different methods or procedures to study 
the effects of reionization history, are respectively referred to as "analytical approximation" , 
"semi-analytical method" and "hydrodynamic simulation". The first is discussed in § ^, 
where a few approximations are made to derive the equation of state in closed form; the 



second is explained in § which evolves the density analytically but requires numerical 
integration of the thermal and chemical evolution equations, though not a full scale hydrody- 
namic simulation and the third is what its name suggests. While the last method is perhaps 
the most accurate (provided the numerical resolution is sufficient), the second, which is the 
focus of this paper, allows us to study the temperature-density relation for a large number 
of reionization histories in an efficient and accurate manner, albeit restricted to the low 
density regime; the first allows us to give useful approximate quantitative expressions that 
embody the correct qualitative trends. Standard symbols are used for cosmological parame- 
ters: a = {l + z)~^ where 1 is the cosmological redshift factor, H for the Hubble constant 
as a function of z, Hq for the Hubble constant today, h for i/o/100kms~^ Mpc^^, for 
the density parameter today, with the subscript b to denote its baryon portion and all its 
matter content, and A for the contribution from cosmological constant. We use the symbol 
h (in distinct from h) to denote the Planck constant in a few places where it arises and /cb 
to denote the Boltzman constant. 



2 THE SEMI-ANALYTICAL METHOD 
2.1 The Equations 

In a nutshell, the method we propose is very simple: assume that the density evolution 
of a given fluid element is governed by the Zel'dovich approximation; allow the element to 
cool/heat and its constituents to change according to standard thermal and chemical reaction 
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rates; repeat the same procedure for a whole set of fluid elements with initial parameters 
drawn from a distribution implied by whatever cosmological initial conditions one prefers (in 
this paper, we study only Gaussian random initial conditions). The resulting temperature- 
density relation can be obtained by making a scatter-plot of temperature versus density for 
all the fluid elements at a given time of interest. 

The density evolution of a fluid element, according to the Zel'dovich approximation, is 
given by: 

1 + 5^ der^[5ij + D+{t)xPij\ , (1) 

where 5 is the mass overdensity, 5ij is the Kronecker delta and -D+(t) is the linear growth 
factor (Peebles 1980), the time dependence of which is completely specified by cosmology 
(e.g. it equals the Hubble scale factor a for a universe at critical matter density). The 
normalization is chosen such that = 1 today. 

The tensor ■0^- is a 3 x 3 symmetric matrix determined by initial conditions (i.e. its 
components vary from one fluid element to another) . One can always choose a basis in which 
the matrix is diagonal. So, the problem of initial conditions reduces to finding the probability 
distribution of its three eigenvalues. This problem was solved by Doroshkevich (1970) for 
Gaussian random initial density field, (see also Bartelmann and Schneider (1992) for an 
elegant alternative derivation) which we also assume in this paper. For recent discussions of 
how to generate realizations of matrix , see Bertschinger and Jain (1994) and Reisenegger 
and Miralda-Escude (1995). It can be shown that the probability distribution of the matrix 
ipij is uniquely specified by the dispersion of its trace. In the limit where -D+(t) is small, 
it can be seen that the trace multiplied by — D+(t) is equal to 5 i.e. the linear overdensity. 
Therefore, the average of the trace vanishes and its rms (root mean squared) value is equal 
to the rms linear density fluctuation today a. At a given redshift and for a flxed number 
of fluid elements, a given rms a value would imply a probable range of 5's in which most 
elements will end up. For the range we are interested in {5 between about —0.9 and 5, at 
z from 2 to 4), we find that setting the rms value a to 2 is adequate. For a fixed number 
of fiuid elements, choosing a higher value for a has the effect of shifting the probable range 
of densities upward and also increasing the scatter in a plot of temperature versus density. 
Alternatively, for a fixed rms value cr, one can reproduce this latter effect and increase the 
number of ffuid elements that fall in the desired density range by increasing the total number 
of ffuid elements in a reahzation. We find that for our densities and redshifts of interest. 
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using a different a does not change the overall temperature-density relation significantly. 
But it should be kept in mind the scatter in temperature-density relation is limited, in part, 
by the finite number of elements in the realizations. 

Temperature obeys the following equation of evolution: 

— = -2HT + _ ^ dE^X. _J_dQ 

dt ?>{l + 5)dt J2iXi dt SkBUb dt ' 

where d/ dt is the Lagrangian derivative following each fluid element, Ub is the proper number 
density of all gas particles (i.e. everything except noninteracting dark matter) and T is 
the temperature which depends on both space and time. The symbol Xi is defined by 
rij = (1 + 6)Xi pb/rrip, where is the proper number density of the specie i, pb is the mean 
mass density of baryons at the time of interest and rrip is the mass of the proton and 6 is the 
mass overdensity. For example the neutral fraction of hydrogen, Xui (to be distinguished 
from Xhi), is then Xhi/ (Xhi+Xhh). Note that Xi is a function of space and time in general. 

The first two terms on the right hand side take care of adiabatic cooling or heating. The 
third accounts for the change of internal energy per particle due to the change in the number 
of particles. The last term dQ/dt is the heat gain (or negative hear loss) per unit volume by 
the gas particles from the surrounding radiation field. The heating and cooling rates due to 
photoionization, recombination and compton scattering are summarized in the Appendix. 

Equation (0) has to be supplemented by one that determines the abundance of each 
particle type, which takes the form: 



—7— — —XiTi + ^ XjXkRjk 



(3) 



where Fj is the photoionization rate of the specie i and Rjk is the recombination rate (in 
units of volume per unit time) of the species j and k to give i. The photoionization rate is 
given by 

where Jp is the specific intensity of the ionizing radiation as a function of frequency Ui is 
the frequency above which a photon can ionize the specie i and ai is the cross-section for 
this process. The cross-sections for different species as well as the recombination rates are 
given in the Appendix. 

What remains to be specified is the time evolution of the specific intensity of the ionizing 
radiation as a function of frequency, which determines the photoionization (Fj) and photoion- 
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Figure 1. A a scatter-plot (black dots) of temperature versus density of 40000 cells randomly drawn from a 64^ CDM+A 
hydrodynamic at z = 4 using the SLH-P''M algorithm (see Gnedin 1995 and Gnedin & Bertschinger 1996) with = 0.35, 
Qa = 0.65, f2f, = 0.055 and h = 0.7 {left panel). The evolution of the ionizing background for the simulation is shown in shown 
in Fig. jij The right panel shows the temperature-density re lati on obtained for exactly the same cosmological and radiation 
parameters using the semi-analytical method described in § 2.1, applied to 2000 fluid elements. The larger scatter in Fig. ^ 
compared to |l^ is in part due to the smaller number of elements in the latter. The solid line plotted in both figures correspond 
to the analytical approximation expressed in equations (p^, (jl^ and (^^ for a sudden reionization model where the epoch of 
rcionization is 2 = 9.6. 

ization heating (dQ/dt) rates. We will examine the effects of a wide range of amplitudes, 
spectra and time evolution of the ionizing radiation in this paper. 

Equations (|l]), (H), and (|^) completely determine the thermal and chemical evolution 
of a fluid element. They are numerically integrated for each fluid element with the initial 
condition for 6 determined as discussed before. The thermal and chemical initial conditions 
are chosen as follows. The gas temperature T is equal to the cosmic microwave background 
temperature at 1 + z = 100(f2{,/i^/0.0125)^/^ (maintained by Compton scattering; Peebles 
1993) and evolves adiabatically after that until the universe is reionized by the UV back- 
ground. Abundances are assumed to be primordial, which is consistent with observations so 
far in the low density intergalactic medium. All species are neutral to high accuracy until 
reionization occurs. One can integrate equations (0) and (|^) forward starting from any time 
between z = 100 and the beginning of reionization. 



2.2 Comparison with Hydrodynamic Simulation 

We show in Fig. |l] scatter plots of temperature versus density at z = 4 obtained using the 
semi-analytical method outlined above and using a full hydrodynamic simulation respec- 
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Figure 2. Evolution of Jhi (solid line), Jhoi {dashed line) and Jhgii {dotted line) as a function of redshift, used in the 
computations shown in Fig. Q The J^'s are defined in equation (^. 

tively. The same uniform ionization field has been used for both of them, the evolution of 
which is shown in Fig. ^.Q The radiation field is computed in a hydrodynamic simulation 
that includes the effects of star formation (see Gnedin and Ostriker 1996). In Fig. |^2, we 
quantify the level of the ionizing radiation by the photoionization rates it implies i.e. 



J, = 10-21 eresHz-^-^cm-^ster-M"^ 



(5) 



4:71 aidu/ V 

where the notations are the same as those defined in equation (|). For instance, Jhi is 
related to the photoionization rate of neutral hydrogen Fhi by Fhi ~ 4 x IQ-^^jj^jg-i (A 
more accurate proportionality constant can be inferred from the cross-section (Thi given in 
our Appendix). The three different quantities Jhi, <^Hei and Jhgh shown in Fig. || do not 
completely specify J^, but they give an indication of how J^, changes with time at different 
frequencies. The full spectrum of Jy is of course known and used in producing the results 



See §1 



for a discussion of the effect of a non-uniform radiation field. 
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shown in Fig. |I] but we do not show it here. (In practice, one actually needs only three 
more pieces of information regarding for the thermal evolution of a gas of primordial 
composition: the photoionization heating rates of Hi, He I and Hell.) 

Our simple semi-analytical method gives a temperature-density relation which is well- 
approximated by a power law for 5^5. It gives the correct mean behavior of gas elements 
in the full hydrodynamic simulation in the low density regime. The larger scatter in Fig. 
|l]a compared to lb is in part due to the smaller number of fluid elements in the latter. At 
higher densities, one begins to see the effect of shocks: a wide scatter of temperature at a fixed 
density. The calculation shown in Fig. |pD is not reliable in this regime. This is because by 
assuming that density evolves as prescribed by the Zel'dovich approximation, which is only 
a good approximation for pressureless (dark) matter or baryons at large scales, one misses 
the effects of gas pressure and shocks as the fluid element is compressed to sufficiently high 
densities. 

We have made a number of similar comparisons between the temperature-density rela- 
tions obtained from hydrodynamic simulations versus using our semi-analytical method, for 
a number of different Jj^'s as a function of time. The semi-analytical method consistently 
gives the correct mean behavior of gas elements of low density. Keeping in mind the intrinsic 
scatter such as that seen in Fig. |l]a, we can make use of our simple semi-analytical method 
to efficiently study the mean temperature-density relation at the low density regime for a 
large number of reionization scenarios, which is the subject of the next two sections. 



3 SUDDEN REIONIZATION MODELS 

3.1 Variation with the Epoch of Reionization 

In Fig. Ha,b,c and d, we show the temperature-density relation for four different sudden 
reionization models, where the reionization epoch is systematically varied. All of them have 
no ionizing background until the specified epoch and then Jm (eq. 0) is taken to be 0.5 

from thereon i.e. 

r Jion for a > Oreion, and 
Jh, = (6) 

10 for a < areion, 

where areion = (1 + -Zrcion)"^ IS the Hubble scale factor when the radiation is turned on and 
the parameter Jion is chosen to be 0.5 (we will study the effects of varying Jjon later on). 
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Figure 3. The temperature-density relation for 4 different sudden-reionization models: sudden reionization (see eq. joj) at 
2 = 5(a),2 = 7(6),2; = 10 (c), and z = 19 (d) . Fo r each reionization model, the black dots shown are the results of calculations 
using the semi-analytical method outlined in § [2.l| for 2000 elements, shown at three different instants; z being 4, 3 and 2 from 
top to bottom. The cosmological parameters are h = 0.5, flo = 1 Eind = 0.0125 (with primordial abundances for hydrogen 

and helium). The ionizing background is specified by its amplitude Jjji = 0.5 (eq. [B|) and spectrum obeying equation M) 
with / = 0.01. The solid, dotted, and dashed lineSjfrom top to bottom for each reionization model, represent the analytical 
expressions for the equation of state in equations (|l^), (^^, and (^2[), using the corresponding cosmological parameters and 
reionization-epoch for each model. The only exce ptio n is panel (d) where the lines shown are exactly the same as those in the 
panel (c), i.e. setting aroion = l/H in equations ( |l9[ ) and (^^ (see explanation in the text). 

It is assumed that Ji, oc z/^^, except that it suffers a factor of / = 0.01 reduction beyond 
the frequency corresponding to the ionization of He II: 

Jo (^) for u < z/Heii, and 

f'^^ (^) > '^Hel, 

where Jq is chosen so that Jhi (eq. |^) has the desired value, uuen is the frequency corre- 
sponding to the ionization of Hell and / (now chosen to be 0.01) is a parameter specifying 
by how much the ionizing radiation is to be further diminished beyond z/hgh. 



J. 



(7) 
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Ignoring for the time being the sohd, dotted and dashed hnes in Fig. |^ which will be 
discussed in the next section, the first trend to notice is that the earlier the epoch of reion- 
ization, the lower the overall temperature at a given redshift, except that for sufficiently 
early reionization (-Zreion > 10) , the temperature-density relation approaches an asymptote. 
The second, somewhat striking, result is how well the temperature-density relations are ap- 
proximated by power-law equations of state, where the slope steepens as universe reionizes 
earlier. The third is that when the reionization epoch is pushed sufficiently early, the exact 
time when it occurs makes very little difference to the equation of state: hence the very 
similar temperature-density relations for reionization starting at 2; = 10 and z = 19 in Fig. 
|c and |d. 

To understand these trends, we perform the following approximate analytical calculation. 
3.1.1 Analytical Calculation 

First, let us determine the reionization temperature. In the limit when a substantial amount 
of radiation is suddenly turned on, one can approximate equations (0) and by: 
dT T dY.,X 2 Xhi /■-. . . , .dv 



dt J2iXi dt 3kBj:iXiJm, 'hz/ 

dXn, V Ta T 

—7— = -Xm / ^T^JuCF-Rir. , 

at JuHi hz/ 

where we have used equation and have taken the crude approximation that the only two 
processes which significantly affect the temperature under the short timescale of reionization 
(reciprocal of photoionization rate) are photoionization heating of hydrogen: 
dQ f°° dp 

-7- ^ nui 47rJi,crHi(hi/ - hz/Hi)— (9) 
dt Ju-iii hz/ 

and the change in the total number of gas particles. Note that during this brief period of 
radiation turn-on, the abundances of various species are far from equilibrium. 

To integrate the above equations, we adopt the further approximation that the universe 
has only hydrogen, i.e. the term ^iXi = 2 — Xhi. Putting T = right before reionization 
because the temperature right after is much higher, we obtain the large time limit of T: 

I 47rJ^(7Hi(hz/ - hz/Hi)— / / AirJ^am — 
i/Hi hz/ Juji, hz/ 

where Ej is kept constant after the radiation field is turned on. This is the temperature that 
the intergalactic medium reaches asymptotically before recombination processes become 
important and halts the exponential decrease of Xhi, setting its value to that implied by 



T ■ = -^Et E 

-'- reion „ 7 ■'-'J 5 
3fcB 



(10) 
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ionization equilibrium thereafter (hydrogen is now highly ionized). Making use of J^, oc u'^, 
it can be shown that Treion ~ 20000 K. This value is not significantly affected by the slope 
of Ju assuming it is between —1 and —1.5 at frequencies immediately above z/hi, nor by the 
behavior of J^, at still higher frequencies. 

Therefore, the first conclusion to draw is that for sudden reionization models, the tem- 
perature of reionization is largely independent of Jhi, density, the epoch of reionization or 
even cosmology. We checked this for a few fluid elements in the computations shown in Fig. 
^ and found that it is indeed the case. One possible exception, which we will discuss in § |3.3| , 
is that changing the spectrum of can affect the helium abundance, which in turn has a 
non-negligible effect on Treion- For sufficiently high redshifts, flreion ~ 1/15, we also find Treion 
to be a little less than the value given above, because of non-negligible amount of Compton 
cooling. 

Next, we consider the evolution of T from the epoch of reionization to 2; = 2, 3 and 4. 
Judging from the result in Fig. ^ that for sufficiently early reionization, the temperature- 
density relation approaches an asymptote, it seems reasonable to assume that what matters 
is the dominant thermal and chemical reaction rates at the later redshifts. In this regime, the 
dominant heating mechanism is again the photoionization of neutral hydrogen, aside from 
adiabatic heating/cooling (recall that we assume a universe filled with only hydrogen). For 
sufficiently low temperature, recombination (of proton and electron) cooling is subdominant 
compared to photoionization heating (at temperature of the order of 10^ K, recombination 
cooling rate is smaller but actually not negligible compared to photoionization heating; it is 
again a crude approximation we adopt to make the problem tractable analytically) . Putting 
all these into equation (|^), the thermal evolution is then approximated by: 

^ = -2HT + + 2 ^ 

dt 3(1 + 5) dt 3kB j:iXi rrip^ ^ ) J \ ) 

where we have ignored the term in equation (|^) due to the change in the number of species, 
used equations (^, (^ and (|TUp and assumed ionization equilibrium (setting the right hand 

side of eq. to zero), with R approximately given by 
T 

i? ^ 4 X 10~^^(Y^)~°-^cm^s"^ (12) 

A more accurate expression for R is given in the Appendix and used for our semi-analytical 
computation but this approximate form suffices for the present analytical calculation. More- 
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over, assuming again there is only hydrogen in the universe, we will approximate Xnn^e/ Hi Xi 
by 1/2 (recall that hydrogen is now highly ionized). 

Even with this highly simplified equation, the problem of predicting the temperature- 
density relation is still untractable analytically, if one uses the density evolution as expressed 
in equation @). One possible approximation is to assume the collapse/expansion occurs 
only along one direction i.e. 1 + 5 = 1/(1 — -D+A) where A is a constant specified by initial 
conditions. Equation ([Tl|) can then be solved analytically but it turns out the resulting 
equations of state differ qualitatively from the ones shown in Fig. ^ in the sense that at high 
and low densities, the predicted temperature tends to be higher than what it actually is. This 
is because in three dimensions, most fluid elements that reach higher densities do so by having 
one axis collapsing while the two others expand. By ignoring the expansion in the other two 
axes (as in a one- dimensional calculation) and thereby missing extra adiabatic cooling, the 
temperature is systematically over-estimated. Moreover, most elements that reach the lower 
densities also have at least two axes expanding and so the one- dimensional calculation again 
underestimates the amount of cooling. The combination of both effects conspire to give an 
equation of state that resembles a power-law for genuine three-dimensional fluid elements. 

The next best question to ask is then: given that the outcome is a power-law equation 
of state, can we at least predict its amplitude and slope? To do so, we turn to linear theory. 
We expand the logarithm of T to flrst order in 5: 

InT = InTo + (7 - 1)5 = InTo + (7 - l)[ln(l + 6)] , (13) 

where Tq is simply the temperature of a fluid element that remains at the cosmic mean 
density and 7 — 1 multiplied by 6 gives the fluctuation about the mean (the notation 7 — 1 



is chosen to agree with the common notation as in eq. [14]). The second equality holds 
in the small 6 limit. The power-law slope relating temperature and density for small 5 is 
then given by dlnT/dln{l + 6) =7 — 1. Knowing that a single power-law equation of state 
approximately holds even for 6 not very small, 7 — 1 is then taken to be the correct power-law 
slope throughout our range of interest i.e. 

T = To(l + 6y-^ (14) 

Putting 6 = into equation (|1T]) and making use of the relation H = HqQo^^'^{1 + z)^/^ 
where Qq is the ratio of matter density to critical density today (the expression is exact for 
f2o = 1 and approximately true at sufficiently high redshift; see Peebles 1993 Chapter 5), it 
is possible to show that 
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(15) 



B 



PclOOm nmO.7 ^6^^ 



T' PT' 

~ rpion -it-t f 



10 K°-^T, 



^ (16) 
_/iv^0.0125_ ^ ' 

where pcioo is the critical density today if = lOOkms"^ Mpc~^, Treion is given by equation 

(p!0|) and R is the recombination rate give in equation ([12|) (note that the combination 

RT^'^ is independent of temperature) and Q^, ^o, h and Hq are cosmological parameters as 

explained in the Introduction. Note also temperature is measured in Kelvin and B has the 

dimension of K^'^, hence the factor of K° '' in the last expression. 



Integrating equation (|15]), we obtain: 
-2 1.7 



rp\.7 / "'reion rp 

V 2 reion 



,1.7 



_3 

H a 2 

1.9 



1 - 



1.9 



B. 



(17) 

reion, the CpOch of 



The initial condition is chosen such that T = Trcion (eq. |]T0[) at a = a 
reionization. 

The equation of evolution for 7 — 1 can be derived from equation ( [TT| ) by expanding T 
as in equation (|13| ) and collecting terms that are first order in 5: 



d{l - 1) 
da 



+ (7 - 1) 



+ 1.7 



d Ina^To 
da 



dhia^Tn 



+ 3 



(18) 



da 6 a 

In deriving the above, we have made use of the fact that the linear growth factor D+, which 
controls the growth of 5 for small 5 {5 (x D^) obeys the relation d\nDj^/d\\\.a = f^o"'^) which 
is exact for critical matter density and approximate otherwise (Peebles 1980). 

Equation (|l^) is a linear ordinary differential equation for 7, which is straightforward, 
if somewhat tedious, to solve. Recalling that Trcion is independent of 6, we use the initial 



condition 7 • 
1 

7 - 1 = — 
' 1.7 



at a 



T 



"reion 
\ 1-7 



and obtain: 



^rcion-^ rcion 



+ 



where C is defined by 



1.9 



,1.9 



1.9 + 1^0 



T 



0.6 



a.9+no" 



2 

3 ~ 
^1.9 



1 

Y.7\ 



[a'Tor'VC, 



0.6 



1 - 



r '^reion xf^o" '^ 

^ a ' 



(19) 



(20) 



\1.7 



''reion reion 



no 



0.6 



,a 



reion ^f^o" 

a 



where B and Tq are defined in equations (|T6|) and (0). The slope 7 — 1 = dlnT / dhi{l + S) 
exhibits the correct qualitative behavior found using our semi-analytical computation shown 
in Fig. 1^. We plot in Fig. ^ the evolution of 7 — 1 as given by equation ([T^) for a number 
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0.6 



0.4 



0.2 




1 + z 

Figure 4. The evolution of 7 — 1 as given by equation ( p^ ) for aroion = 1/6 {solid line), a^, 
iroion = 1/11 {dashed line). Here Troion = 20000 K and flo = 1 in all cases. 



1/8 {dotted line) and 



of different ctreion- An interesting consequence of the above expression for the slope is a hmit 
as to how steep the temperature-density relation can be (or an upper limit to 7 — 1). In the 
limit where reionization occurs very early: 



7i 



- 1 



1 0.078^0 

+ 



1.7 ■ i.g + no^*^- ^^^^ 

Hence, the maximum possible value for 7 — 1 is 0.62 (assuming Qq ^ !)• This has interesting 
implications for the slope of the column density distribution of the Lyman alpha forest, which 
we will discuss in § |^. 

We show in Fig. |^ a comparison of the above expressions (eq. [1^ and |19[) with the 
temperature-density relation obtained by solving numerically the rate equations as outlined 
in Sec prT| . The evolution of the slope is quite accurately captured by our analytical approx- 
imation in equation (plQl). However, the amplitude Tq predicted by equation ([17|) is off: the 
overall predicted temperature is lower than what it actually is and the the predicted spread 
in To from z = 2 to 2; = 4 is also larger than what the semi-analytical calculation shows (Fig. 
^). This is mainly because helium is not taken into account at all in the above treatment. 
By ignoring heating due to photoionization of helium, we systematically underestimates Tq 
in our analytical calculation. The time evolution for Tq is also not quite right for the same 
reason. We therefore introduce the following modification of equation (|l^) that gives an 
excellent fit to the results of the semi-analytical calculations shown in Fig. |^: 
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Figure 5. Temperature-density relation for two sudden-reionization models: comparison of results of semi-analytical calcula- 
tions (black dots) and predictions based on analytical approximations expressed in equations (^), {[^, and Q {solid, dotted 
and dashed lines). The model parameters are described in the caption of Fig. ^. The black dots in the left panel and right panel 
are exactly the same points as those in Fig. ^ and Fig. ^ respectively , i.e. sudden reionization at ^ = 6 and z = 19. From 
top to bottom for each figure, the results are shown at z = 4, 3, and 2 respectively for the black dots and also for the lines. 
Notice how the predicted slopes match quite well while the amplitudes are off. A modified version of equation ( p^ ) (eq. p2| ) 
gives much better fit to the amplitudes and is used instead in the rest of this paper. 



rpl.7 / '^reion rji 

\ o reion 



M7 1-7 



3 I a-0.25 
' 2 a 



1.9' 



B 



T 

J- rf 



25000 K 



(22) 



where we have modified the power index of a in the second term on the right hand side and 
used Treion = 25000, in place of Tj-eion given in equation (|TD|). In general, a different best-fit 
7"reion would have to be chosen for different abundances of the helium species, when one 
varies the spectrum of the ionizing radiation, for instance. We will return to this point later. 

On the other hand, the slope of the equation of state 7 — 1 is quite insensitive to the 
exact value of Treion, and using the original estimate (eq. |jlO|) in equation ([T9| ) is sufficiently 
accurate. 

The analytical predictions for the equation of state shown in Fig. ^ are based on this 
modified version of Tq (eq. |2^) together with equations (|13D and (0). Note that, for the 
case in which the universe reionizes at z = 19, the temperature-density relation is very close 
to that of the case in which the universe reionizes at z = 10. This turns out to be a general 
trend, that there is an asymptotic temperature-density relation when the universe reionizes 
sufficiently early. To see how this asymptotic state is approached, we show in Fig. |1^ the 
evolution in temperature of one particular fiuid element for different reionization histories. 
Notice how the solid line (sudden reionization aX z = 10) and the dotted line (sudden 
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Figure 6. The temperature-density relations for Qi,h^ = 0.025 (upper black dots) and Qi,h^ = 0.0125 (lower black dots) 
at z = 3 for sudden reionization model with radiation-turn-on at 2 = 7. The rest of the cosmological parameters are exactly 
the same as those given in Fig. jsj The dotted and solid lines are the predictions according to equations (p^, ( p^ ) and (p^), 
respectively for the two models. 

reionization at z = 19) converges. The other hnes in the figure corresponds to reionization 
models that will be discussed later. We also find that fiuid elements at other densities also 
exhibit the same kind of convergence towards an asymptotic state. 

The expressions in equations (|1^), (p^, and (|22|) are derived by ignoring processes like 
Compton cooling which are important at higher redshifts, if the universe is already reionized 
by then. Hence, the analytical expressions are not quite as accurate for areion < 1/11- How- 
ever, we can make use of the knowledge of asymptotic behavior for early reionization models 
and use simply the areion = 1/11 result of equations (|13]), (pJ]), and p^ ) for any models in 
which reionization occurs earlier than z = 10. This is what we have done in making Fig. ^d. 

We now turn to two other interesting variations of the sudden-reionization models. 



3.2 Variation with Cosmology 

We investigate the temperature-density relations for different cosmological parameters, in- 
cluding the baryon density, the Hubble constant and the fraction of critical density in mat- 
ter. One example is shown in Fig. ^, where Qf,h'^ has been increased by a factor of two from 
the value in Fig. |^. Once again, we find that a power-law equation of state gives an ade- 
quate description of the temperature-density relation. Moreover, the analytical expressions 
in equations (0), ([T9|), and (^2]) work quite well in describing the change in amplitude and 
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slope. As predicted by equations (122|) and (p^Bf ), the overall temperature increases as Qbh'^ 
is raised because the rate of photoionization heating is increased. (Recombination cooling 
rate is actually also increased but photoionization heating is the dominant entropy changing 
process.) 

We have also tried other variations such as lowering Q^. A realistic structure formation 
model such as the cosmological constant dominated Cold Dark Matter model with Qq = 0.35 
and h = 0.7 (Kofman, Gnedin, & Bahcall 1993; Ostriker & Steinhardt 1995) does not give 
a very different temperature-density relation from the ones shown in Fig. |^ (using the same 
Qbh"^). This can be understood from equations (^) and (0), where it is clear the overall 
temperature to the 1/1.7-th power is dependent upon the combination h^/^l^ which only 
varies from 0.5 for the model in Fig. ^to about 0.4 in the cosmological constant model. 

In general, we find that equations (|13]), ( |T9| ) and ( p2D give an agreement to the actual 
temperature-density relation to within 20% for a variety of cosmological models we have 
tested. A large part of this error has to do with the fact that changing Qq and h changes the 
amount of time helium photoionization heating or recombination cooling has to act to make 
a difference to the thermal evolution. The effect of helium is what we will discuss next. 

3.3 Variation with Amplitude and Spectrum of Radiation 

The last set of tests we perform for the class of sudden reionization models is to vary the 
spectrum and amplitude of the ionization radiation. We find that varying the amplitude 
of the ionization radiation does not affect the temperature-density relation significantly, 
as predicted by equations (|I9|) and (p2D, as long as the flux level is sufficient to maintain 
the universe highly ionized. This is shown in Fig. ^ and |^, which corresponds to choosing 
Jion = 0.1 and Jion = 2.0 respectively in equation (^). A comparison of the two figures reveals 
a slight decrease in temperature when Jhi is raised, which is a result of more opportunities 
for recombination cooling, but the effect is very small and becomes insignificant at lower 
redshifts {z = 2). 

The main effect of varying the spectrum (when fixing Jhi defined in eq. [H) is that one 
changes the relative abundance of the different species of helium, which could have a non- 
negligible effect on the temperature due to recombination cooling of the ionized species and 
photoionization heating of neutral helium. We test this effect by varying the Hell cut-off 
factor / in equation (|^). Two examples are shown in Fig. ^ where we replot the lines in Fig. 
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Figure 7. The effect of varying the amphtude of Ji, on the temperature-density relation. Black points in the left panel (Fig. 
^) are the result of semi-analytical calculation of sudden reionization at a^eion = 1/8 with Ji^n = 0.1 (eq. [Q) while those in 
the right panel (Fig. ^) has Ji^n = 2.0 starting from the same arcion- The spectrum is chosen to be the same with / = 0.01 
(eq. M) for both panels. The outputs from top to bottom are at ^ = 4, 3 and 2. The same solid, dotted, and dashed lines are 
drawn for both panels to facilitate comparison. They are also exactly the same as the lines in Fig. ^b, based on the analytical 
expressions in equations (p^), (p^, and (p2[), using arcion = 1/8- The cosmological parameters in Fig. ^ are exactly the same 
as those in Fig. 
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Figure 8. The effect of varying radiation spectrum on the temperature-density relation. Black dots are the results obtained 
by the semi-analytical calculation explained in § ^.l| . The reionization and cosmological model here is, except for the spectrum, 
exactly the same as that in Fig. |3^ where the universe reionizes suddenly at z = 6. The solid, dotted and dashed lines of Fig. ^ 
are reproduced here for comparison. Above and below the solid line at the top are two sets of points at 2 = 4 which correspond 
to choosing / = 1.0 (less Hell) and / = 10"^ (more Hell) respectively in equation (^. The same is also true for the two sets of 
points surrounding the dotted line, except that they are at z = 3. Results are not shown for z = 2 because the two different /'s 
give significantly overlapping results that are well-represented by the dashed line. The amplitude Jm = 0.5 (eq. ^) is adopted 
throughout. 
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^ (which match the temperature-density relations for / = 0.01 quite well) together with 
the result of two semi-analytical calculations for / = 1.00 and / = 10^^ respectively. By 
decreasing / or the amount of radiation energetic enough to ionize He ll, one increases the 
abundance of He ll offering more opportunities for recombination cooling, thereby lowering 
the overall temperature. The reverse happens when / is raised. The change in temperature 
due to this effect could be as large as 20%, with less change as one goes to lower redshift, 
for a fixed reionization history. One can compensate the error in the analytical estimate 
of To (which ignores the effect of helium) by lowering/raising Treion in equation (^2]) as 
one decreases/increases / (For instance, choosing Treion = 28000 for the / = 1 case and 
^reion = 22000 for the / = 10~^ works reasonably well). 

The power-law slope, on the other hand, does not vary significantly with the spectrum of 
the ionizing background. It can be seen from Fig. ^ that there is a slight decrease of slope for 
the corresponding redshifts as one lowers /. This is due to the fact that Hell recombination 
cooling has a stronger effect at higher density. The analytical prediction for the slope in 
equation (|T9|) works best when the effects of helium is negligible, as can be seen from the 
better match of slopes for the / = 1 case. Nevertheless, on the whole, we find that the 
prediction of equation (|T^ for the slope is quite robust. 

To summarize the results of this section, we find that the temperature-density relation 
at 2 of 2 to 4 for sudden-reionization models is well described by a power-law equation 
of state, with a scatter that depends on the total number of fluid elements. The overall 
amplitude decreases and the slope steepens as the epoch of reionization is pushed earlier; 
but the equation of state rapidly approaches an asymptotic limit for reionization occurring 
at redshifts beyond 10. Quantitative estimates for the slope and amplitude are given in 
equations (|19D and (|22|). While the analytical prediction for the slope is relatively robust 
for various cosmologies and spectra, the prediction for the amplitude agrees with the results 
of the semi-analytical calculations (outlined in § |2]2) to within 20% for reasonable range of 
cosmological and ionization parameters. Most of the error comes from the neglect of helium 
in obtaining the analytical estimate, which can be compensated if one allows Tj-oion to change 
in equation ( P^ according to the abundance of the helium species (Fig. 
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4 REIONIZATION PRECEDED BY REHEATING 

The idealization of sudden reionization in the last section allows us to understand the various 
trends of the temperature-density relation analytically. It is quite possible that the ionizing 
radiation is turned on more gradually in nature. In particular, we focus on a class of models 
in which radiation is turned on at some early time, with an intensity large enough to heat 
up the intergalactic medium but small enough so that the neutral fractions remain high 
(large recombination rates) though slowly decreasing. This is motivated by hydrodynamic 
simulations (see Gnedin & Ostriker 1996) in which collapsed regions of high density (could 
be stars/quasars) are assumed to give out radiation. As the formation rate of these regions 
increases with time and the neutral fraction drops steadily, the ionizing intensity climbs until 
at some point, the universe becomes highly ionized (neutral fraction of hydrogen less than 
10~^), thus completing the reionization process when the radiation levels off and species 
abundances are simply given by ionization equilibrium. The simulation shown in Fig. ^a 
and ^ is one such example. 

To investigate the implications of these reionization models systematically, we adopt the 
following parametrization for the evolution of the ionizing radiation: 

{Jion for 
-''heat (^) for ttion > a> ttheat, (23) 
for ttheat > a, 

where Jhi is defined in equation (||), ahcat and aion denote the size of the Hubble scale factor 
at the onset of reheating and reionization. The slope a'' in the second row of equation (^31) 
comes from fitting Jm{(i) in Fig. ^ The spectrum of is assumed to have the same form 
as in equation (0) with / = 0.01 and Jq chosen to match whatever Jhi specified above. 

Before discussing the results for these models, let us point out that in Fig. ^ we have 
already shown an example of the temperature-density relation of a reionization-preceded-by- 
reheating model where J varies with time as in Fig. ^. Notice how the temperature-density 
relation can be approximated by an equation of state resulting from a sudden reionization 
model (the solid line in Fig. |I]). The same turns out to be true for all other reionization 
models we have tested, with the parametrization given in equation (PBf). 

Two interesting examples are shown in Fig. ^. Both have the same epochs of reheating and 
reionization (oheat and Oion)- The only difference is the amount of reheating that occurs: Jheat 
being 0.001 versus 0.1 in equation (p3|). They can both be approximately fitted with sudden 
reionization models (no preceding reheating period) where the epoch of sudden reionization 
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Figure 9. Black dots in Fig. ^a on the left show the temperature-density relation for a reionization model in which ajieat = 1/20, 
lion = 1/6, Jhcat = 0.001 and Jjon = 0.5 (see eq. ]23|). The three sets of points are from top to bottom at z = 4, 3 and 2. The 
solid, dotted and dashed lines are analytical approximations of sudden reionization models expressed in equations (p^), (p^), 
and 023), choosing Qrcion = 1/8.4. Fig. tip on the right, on the other hand, corresponds to the choice of parameters aj^eat = 1/20, 
gjnn =1/6. JhGaj_= 0.1 and Jion = 0.5 in cq. [ pif . The analytical approximations shown have arejon equal to 1/10 in equation 
(|l3[), ( p^ and (p^). For both Fig. ^a and the cosmological model is the same as that in Fig.^ 

is chosen to lie somewhere between Cheat and aion- The larger the amount of reheating, the 
earlier the epoch of reionization is required. It is also interesting to note that in Fig. |^a, the 
fit to the result at z = 4 is not as good as those at 2; = 3 or 2; = 2. It is because it takes some 
time for the universe to settle to its asymptotic state (i.e a state which is independent of 
prior history, whether it be sudden reionization, gradual reionization with preheating, etc). 
The same should also be true for Fig. Pd, except that because the amount of preheating 
is more significant, the universe, in effect, has more time to settle to its asymptotic state, 
hence the better agreement with sudden reionization models for all three redshifts. 



To illustrate this convergence towards an asymptotic state, we show in Fig. |10| the evolu- 
tion of temperature for one particular fluid element for four different reionization histories. 
The tendency to approach the same limiting temperature, for a given density evolution, 
can be clearly seen. Note how the model illustrated with a long-dashed line evolves to a 
temperature a little higher than the others, by z = 2. This is because of its relatively small 
amount of reheating (small Jheat in eq. ||23|) which implies the universe effectively reionizes 
at a latter redshift, giving it less time to settle to its asymptotic state. 

We have tested several other reionization models where all the parameters in equation 
(p3|) are varied systematically. The same conclusion holds for all of them, namely that a sud- 
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Figure 10. Examples of convergence towards an asymptotic state. In the upper panel is shown the evolution of the overdensity 
of a fluid element as a function of redshift. In the lower panel is shown the evolution of temperature for the same fluid element, 
but for a variety of reionization histories. The solid and dotted lines both correspond to sudden reionization models in which 
o-ieiou = 1/11 S'lid <ireion = 1/20 respectively and Ji^n = 0.5 in equation Their reionization and cosmological parameters 
are exactly the same as those in Fig. He and bp respectively. The short dashed and long dashed lines correspond exactly to 
the reionization models shown in Fig. Hb and ffi respectively i.e. reionization-preceded-by-reheating models with aha at = 1/20, 
lion = 1/6, Jion = 0.5 and Jhcat = 0.1 for short dashed line and Jheat = 0.001 for long dashed line (see eq. ]23[). The 
cosmological parameters for these two latter models are the same as those for the solid and dotted lines. 

den reionization model (no preceding reheating) can always be found to fit the temperature- 
density relation of any reionization model with reheating, with the turn-on redshift of the 
former chosen somewhere between the epochs of reheating and reionization of the latter. 

5 DISCUSSION 

Let us summarize what we have learned so far. 

We find that the temperature-density relation of the low density intergalactic medium, 
while having an intrinsic scatter, is well described by a mean equation of state of the form 
T = To(l + 6y'~^. It is shown that any reionization model with gradual radiation-turn-on 
produces equation of state that can be well approximated by one produced by a suitably 
chosen sudden turn-on model (§ 

In general, the mean temperature Tq at 5 = at a given redshift ^ ~ 2 — 4 decreases 
as the epoch of reionization is pushed earlier. It is only weakly dependent on the ampli- 
tude of the ionizing flux Jhi at z ~ 2 — 4 assuming reionization takes place at z^5. The 
spectrum of can introduce about 20% change to Tq. Its main effect is on the abundances 
of helium species, affecting the amount of photoionization heating and recombination cool- 
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ing. Cosmological parameters change Tq in a simple way: for sufficiently early reionization, 
To oc {^Ibh/y/Tl^y^^-'^ . All the dependence mentioned above are expressed in a quantita- 
tive way in equation (p2D, with about 20% accuracy depending on the spectrum of Ju and 
cosmology. 

The slope of the equation of state, 7 — 1, increases with redshift and with earlier reion- 
ization epoch (eq. [jl9l and Fig. ||). Its evolution is only weakly dependent on cosmology. 
There exists a maximum as to how steep 7 — 1 can get, which is about 0.62 (eq. 
Strictly speaking, 7 — 1 can be as small as right at the onset of sudden reionization, but 
for reionization occurring earlier than Zreion ~ 5, it is greater than 0.3 at z = 3; for ^reion > 7 
the minimum 7 is about 0.47. 

What do all these imply for the observable properties of the low column density Lya 
forest? Firstly, for a given density (or 1 + 6) profile and the Hubble constant h, the column 
density it corresponds to is proportional to ^IfT^^'"^ J^] . Uncertainty in the reionization 
history introduces about 50% uncertainty in T (see Fig. which in turn implies about 30% 
uncertainty in column density. This means an additional uncertainty in the amplitude of the 
column density distribution on top of that due to fit and Jhi, even if one is given a known 
cosmological model (Hui et al. 1996). 

Note that for sufficiently early reionization, the dependence Tq oc ^b^^^''^ implies that the 
column density is scaled by ■ ' ' = which is roughly consistent with the result of 
Croft et al. (1996), who found Ql'"^ scaling using hydrodynamic simulations. 

The good news, on the other hand, is that the slope of the column density distribution is 
relatively independent of reionization history. The slope of the column density distribution 
(number of absorption lines per unit column density per unit redshift oc A^nf where N^i is 
the column density) is approximately given by: 

where cxo is the rms mass fluctuation at scales close to the Jeans length. It is close to 1 but 
its precise value is determined by the cosmological model (Hui et al. 1996). Taking the range 
of 7 — 1 from 0.3 to 0.62 we have derived before and Uo ~ 1, the uncertainty in (3 introduced 
by unknown reionization history is of the order of 5%. 

This is a very important conclusion because the slope of column density distribution then 
becomes a powerful discriminant of cosmological models through a^. It is relatively free of 
uncertainties due to reionization history, or the values of fib and Jri, unlike the amplitude 



© 0000 RAS, MNRAS 000, 000-000 



26 Hui and Gnedin 

of the column density distribution. The reader is referred to Hui et al. (1996) for a more 
detailed discussion of the predictions of various cosmological models for p. 

The temperature-density relation also has significant implications for the b parameter 
distribution. There exists a rather sharp lower cut-off of about b = 20kms~^ according to 
high resolution data of the low column density Lya observation (Hu et al. 1996). If one 
naively assumes the observed b is completely a result of thermal broadening b = ^2kT/mp, 
one obtains a lower cut-off temperature of about 24000 K. Judging from Fig. ^, it is clear 
that at 2 ~ 3, which is the relevant redshift for these observations, unless the universe 
reionizes quite a bit later than z = 5, all the underdense regions (and including some 
overdense regions too, depending on the reionization history) have lower temperature that 
this cut-off. This is true even if one allows for a lot of helium photoionization heating or very 
little helium recombination cooling by varying the spectrum (see Fig. There are a few 
possible resolutions: first, it is possible that the low column density objects that we observe 
(column density lower than about 10^^ cm^) actually correspond to higher densities (6) , 
but that would imply a very high Jhi (higher than about 2.0) which is difficult to reconcile 
with observations (see Hui et al. 1996); second, it could be that the universe is not photo- 
reionized after all and some other mechanism is at work; third, the observed b parameter 
for these low column density systems is not the result of thermal broadening alone. We 
consider the last possibility the most likely one. As we have argued in an earlier paper (Hui 
et al. 1996) and was pointed out before by numerous authors (see for instance Cen et al. 
1994; Hernquist et al. 1996; Miralda-Escude et al. 1996; Zhang et al. 1995), the observed 
absorption profiles do not necessarily conform to the classic Voigt profile but can instead 
reflect more the underlying mass density profile and peculiar velocity structure. For instance, 
a structure collapsing perpendicular to the line of sight (thus having density enhancement 
to induce enhanced absorption) can have an expanding velocity component along the line 
of sight, thus creating a profile broader than what naive thermal broadening would predict. 

Finally, let us speculate on how a fluctuating ionizing field might change the results of 
our paper. From Fig. ^, it is clear that a fluctuating field with inhomogeneous spectrum can 
certainly introduce a lot of scatter to the temperature-density relation, by giving different 
photoionization heating and recombination cooling rates even for fluid elements with the 
same density evolution. Moreover, the ionizing field can also turn on at different times 
for different fluid elements because of radiative transfer effects (Zuo 1992a,b). This would 
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introduce additional scatter because as we have seen, the earher the reionization epoch, the 
lower the temperature at a given redshift. 

However, we also learn that the temperature-density relation approaches an asymptotic 
limit if one waits long enough. For sufficiently early reionization, neither the amplitude of 
Jhi nor the reionization history (eg. at what different times radiation turns on at different 
elements) affects the temperature very much, assuming Jhi is large enough to keep the 
universe reionized (see § p.3|) . Therefore, provided the amplitude of Ju at the helium ionizing 
frequencies does not vary significantly with position at later times z ~ 2 — 4, the equation 
of state we calculate for the early reionization models should be valid, even if the amplitude 
of Jhi fluctuates at 2; ~ 2 — 4. This certainly warrants further research, which we will leave 
for future work. 
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APPENDIX A: IONIZATION AND RECOMBINATION RATES 

In this section we compile the ionization and recombination rates used in our calculations. 
We believe that this compilation is sufficiently up-to-date to reflect all recent improvements. 
In some of the equations below the following notation is used: 



where T is the temperature, and Tp'^^ are ionization thresholds for species j = H I, He I, He ll 
expressed in temperature units, 
T^™) = 157807 K, 
T^f^ = 285335 K, 
= 631515 K. 
The symbol ks denotes the Boltzmann constant. 

For chemical evolution, we use the symbols RI, CI and DI to denote hydrogen and 
helium recombination rates, coUisional ionization rates and dielectronic recombination rates. 
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For thermal evolution, we use the symbols RC, CC, DC and EC to denote hydrogen 
and helium recombination cooling, coUisional ionization cooling, dielectronic recombination 
cooling and line excitation cooling rates. 

Superscripts and subscripts are added to the above symbols to distinguish between dif- 
ferent processes within each class. 

Molecular hydrogen and metal cooling are not expected to be important for the low 
density intergalactic medium that we are interested in (see Songaila & Cowie 1996 for a 
discussion of metals detected in the Lya forest). 

Compton cooling rates and the photoionization cross-sections are given at the end. 

case A Hll recombination coefficient: our fit to the data from Ferland et al. (1992); accu- 
rate to 2% from 3 K to 10^ K: 

\ 1.503 

RI±, = 1.269 X lO-^^cm^sec"^ Si 

(1.0 + (Ah,/0.522)0-470)1-923 

case A Hll recombination cooling rate: our fit to the data from Ferland et al. (1992); ac- 
curate to 2% from 3K to 10^ K: 

\ 1.965 

RC4u = 1.778 X lO-^'^ergcm^sec-^K-^T ^ ™ 

(1.0 + (Ah,/0.541)0-502)2-697 

case A Hen recombination coefficient: from from Burgess & Seaton (1960); accurate to 
~10% from 5 X 10^ K to 5 X 10^ K: 

= 3.0 X lO-^^cm^sec-^A^^?^ 

case A Hell recombination coohng rate: from from Burgess & Seaton (1960); accurate to 
~10% from 5 X 10^ K to 5 X 10^ K: 

case A Hem recombination coefficient: our fit to the data from Ferland et al. (1992); ac- 
curate to 2% from 1 K to 10*^ K: 

\ 1.503 

RI±,,, - 2.0 X 1.269 X 10-^2 cm^ sec"^ Sen 

(l.O+(AHen/0.522)0-470)l-923 

case A Hem recombination coefficient: our fit to the data from Ferland et al. (1992); ac- 
curate to 2% from 1 K to 10^ K: 

\ 1.965 

^C'Hem = 8 X 1.778 X lO-^^ergcm^ dimsec-^K'^T ^ ™ 

^ (1.0 + (AHe„/0.541)0-502)2-697 

case B Hll recombination coefficient: our fit to the data from Ferland et al. (1992); accu- 
rate to 0.7% from 1 K to 10^ K: 
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\1.500 

Rli„ = 2.753 X ^ ™ 

Hn (1.0 + (Ah,/2.740)0-407)2-242 

case B Hll recombination cooling rate: our fit to the data from Ferland et al. (1992); ac- 
curate to 2% from 1 K to 10^ K: 

\ 1.970 

RCi,, = 3.435 X 10-^° erg cm^ sec"^ K'^T 



,1.0 + (Ah,/2.250)0-376)3-720 

case B Hell recombination coefficient: from from Burgess & Seaton (1960); accurate to 
-10% from 5 X 10=^ K to 5 X 10^ K: 



= 1.26 X lO^i^cm^sec-U^Jf 



case B Hell recombination cooling rate: from from Burgess & Seaton (1960); accurate to 
~10% from 5 X 10^ K to 5 X 10^ K: 

case B Hem recombination coefficient: our fit to the data from Ferland et al. (1992); ac- 
curate to 2% from 3K to lO'^ K: 

X 1.500 

B _ o n o TCQ in-14 „™3„„„-l ^Hcii 



i?/^^,„ = 2.0 X 2.753 X 10"^^ cm^ sec"^ ™ 

(1.0 + (AHen/2.740)0-407)2-242 

case B Hem recombination coohng rate: our fit to the data from Ferland et al. (1992); 
accurate to 2% from 3K to 10^ K: 

\ 1.970 

RCi^m = 8 X 3.435 X 10-^° erg cm=^ sec"^ K'^T 



1.0 + (AHei,/2.250)0-376)3-720 

Hi collisional ionization coefficient: our fit to the data from Lotz (1967); accurate to 3% 
from 10^ K to 10^ K: 



C7hi = 21.11 cm^sec-^K=^/2r-^/2e-^Hi/2 ^ 



1.089 



:i + (Ah:/0.354)0-874)1-101 

H I collisional ionization cooling rate: derived from the collisional ionization coefficient: 

He I collisional ionization coefficient: our fit to the data from Lotz (1967); accurate to 3% 
fi-om 10^ K to 10^ K: 

, \-1.146 

C/Hei = 32.38 cm=^sec-iK=^/2T-=^/2e-^Hei/2 ^Hc^ 

(1 + (AHei/0.416)0-987)l-056 

He I collisional ionization cooling rate: derived from the collisional ionization coefficient: 
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Table Al. Photoionization cross-section parameters 



Species 


Eo (eV) 


o-Q (cm^) 


P 


Va 






yi 


Hi 


4.298 X 10-^ 


5.475 X 10-^"* 


2.963 


32.88 











He I 


1.361 X 10^ 


9.492 X 10-1^ 


3.188 


1.469 


2.039 


0.4434 


2.136 


Hen 


1.720 


1.369 X lO-i"* 


2.963 


32.88 












Hell collisional ionization coefficient: our fit to the data from Lotz (1967); accurate to 3% 
from 10^ K to 10^ K: 

\ /o \ -1.089 

C/He„ = 19.95 Cm3seC-lK3/2T-3/2e-AHen/2 ^Hcn 

(l + (AHen/0.553)0-735)l-275 

Hell collisional ionization cooling rate: derived from the collisional ionization coefficient: 

CC*Hcn = ^B^hIh'' C*-^Heii 

Hell dielectronic recombination coefficient: from Aldrovandi & Pequignot (1973); accurate 
to ~5% from 3 x 10^ K to 10^ K: 

D/hcii = 1-90 X 10~3cm3sec-^K3/2r-3/2e~0-''^^Heni/2 ^ Q_3g-0.15AHcin/2^ 

He II dielectronic recombination cooling rate: derived from the dielectronic recombination 
coefficient: 

-DChcii = 0.75A;B^Hlm -D/hchi 

Hi line excitation cooling rate: from Black (1981) with the correction from Cen (1992); 
accurate to ~10% from 5 x lO^K to 5 x 10^ K: 

EChi = 7.5 X 10-^^ergcm=^sec-^e"0-''^^Hi/2 / 

& l + (r/105K)V2 

Hell line excitation cooling rate: from Black (1981)with the correction from Cen (1992); 

accurate to ~10% from 5 x 10^ K to 5 x 10^ K: 

EChch = 5.54 X 10-^ erg cm^ sec- j^l^^'^^^ ^-0.75AHe„/2 ^ ^ ^^^j^^ ^^^^^ 

Compton heating/ cooling term in the equation (D) is given by the following expression 
(Peebles 1993), which is exact in the nonrelativistic limit: 

c^gcompton ^ g 35 ^ iq-4i ^^-3 ^-1 ^-i^^h^x,{l + (2.726 K{1 + z) - T) . 



Photoionization cross-sections are taken from Verner et al. (1996) in the following form: 

0.5P— 5.5 



y 

1 + Jy/va 



p 1 
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where x = E/Eq — yo and y = \Jx^ + yf. Quantities (Tq, Eq, y^, P, ya, yo, and yi are 
fitting parameters and are given in the Table 0. They are accurate to within 10% from the 
respective ionization thresholds to 5 keV. 



REFERENCES 

Aldrovandi, S. M. V., & Pequignot, D. 1973, A&A, 25, 137 
Bartelmann, M., & Schneider, P. 1992, A&A, 259, 413 
Bertschinger, E., & Jain, B. 1994, ApJ, 431, 486 
Bi, H. G., Borner, G., & Chu, Y. 1992, A&A, 266, 1 



Bi, H. G., & Davidson, A. F. 1996, ApJ, in press ( ^stro-ph 9611062[ ) 

Black, J. H. 1981, MNRAS, 197, 553 

Burgess, A., & Beaton, M. J. 1960, MNRAS, 121, 471 

Cen, R. 1992, ApJS, 78, 341 

Cen, R. Y., Miralda-Escude, J., Ostriker, J. P., & Ranch, M. R. 1994, ApJ, 437, L9 
Coles, P., Melott, A. L., & Shandarin, S. F. 1993, MNRAS, 260, 765 
Couchman, H. M. P., & Rees, M. J. 1986, MNRAS, 221, 53 



Cristiani, S., D'Odorico, S., D'Odorico, V., Fontana, A., Giallongo, E., & Savaglio., S. 1996, MNRAS, in press (istro-ph 



9610006) 



Croft, R. A. C., Weinberg, D. H., Katz, N., & Hernquist, L. 1996, ApJ, submitted ( |astro-ph 9611053| ) 
Doroshkevich, A. G. 1970, Astrofizika, 6, 581 

Ferland, G. J., Peterson, B. M., Horne, K., Welsh, W. F., & Nahar, S. N. 1992, ApJ, 387, 95 
Giroux, M. L., & Shapiro, P. R. 1996, ApJS, 102, 191 
Gnedin, N. Y. 1995, ApJS, 97, 231 

Gnedin, N. Y., & Bertschinger, E. 1996, ApJ, 470, 115 
Gnedin, N. Y., & Hui, L. 1996, ApJ, 472, L73 



Gnedin, N. Y., & Ostriker, J. P. 1996, ApJ, submitted ( ^stro-ph 9612127| ) 
Gunn, J. E., & Peterson, B. A. 1965, ApJ, 142, 1633 
Haiman, Z., & Loeb, A. 1996, ApJ, submitted ( astro-ph 9611028 ) 
Haardt, F., & Madau, P. 1996, ApJ, 461, 20 

Hernquist, L., Katz, N., Weinberg, D. H., & Miralda-Escude, J. 1996, ApJ, 457, L51 

Hu, E., Kim, T., Cowie, L. L., & Songalia, A. 1995, AJ, 110, 1526 

Hui, L., Gnedin, N. Y., & Zhang, Y. 1996, ApJ, submitted ( astro-ph 9608157 ) 

Ikcuchi, S., & Ostriker, J. P. 1986, ApJ, 301, 522 

Kofman, L. A., Gnedin, N. Y., & Bahcall, N. A. 1993, ApJ, 413, 1 

Lotz, W. 1967, ApJS, 14, 207 

Lu, L., Sargent, W. L. W., Womble, D. S., & Takada-Hidai, M. 1996, ApJ, 472, 509 
Meiksin, A.,& Madau, P. 1993, ApJ, 412, 34 

Miralda-Escude, J., Cen, R., Ostriker, J. P., & Ranch, M. 1996, ApJ, 471, 582 
Miralda-Escude, J., & Ostriker, J. P. 1992, ApJ, 392, 15 
Miralda-Escude, J., & Rees, M. J. 1994, MNRAS, 266, 343 
Ostriker, J. P., & Gnedin, N. Y. 1996, ApJ, 472, L63 
Ostriker, J. P., & Steinhardt, P. J. 1995, Nature, 377, 600 

Peebles, P. J. E. 1980, The Large Scale Structure of the Universe, Princeton: Princeton University Press 



© 0000 RAS, MNRAS 000, 000-000 



32 Hui and Gnedin 



Peebles, P. J. E. 1993, Principles of Physical Cosmology, Princeton: Princeton University Press 
Reisenegger, A., & Miralda-Escude, J. 1995, ApJ, 449, 476 
Sciama, D. W., 1990, ApJ, 364, 549 
Shapiro, P. R. 1995, ASP Conf. Ser., 80 

Shapiro, P. R., Giroux, M. L., & Babul, A. 1994, ApJ, 427, 25 
Songaila, A., & Cowie, L. L. 1996, ApJ, 112, 335 

Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlcv, D. G. 1996, ApJ, 465, 487 

Wadsley, J. W., & Bond, J. R. 1996, preprint ptro-ph 9612"l48 , to appear in Proc. 12th Kingston Conference, Halifax 
Zel'dovich, Ya. B. 1970, A&A, 5, 84 

Zhang, Y., Anninos, P., & Norman, M. L. 1995, ApJ, 453, L57 
Zuo, L. 1992a, MNRAS, 258, 36 
1992b, MNRAS, 258, 45 



© 0000 RAS, MNRAS 000, 000-000 



